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The quantum theory of the Coulomb field has been developed by Staruszkiewicz in the long series 
of papers. This theory explains the universality and quantization of the electric charge observed in 
Nature. Moreover, the efforts have been made to determine the value of the elementary charge from 
its mathematical structure. Nonetheless, no other immediate applications of this theory have been 
■ proposed. We make such an attempt by (i) considering the classical energy operator and defining 

its counterpart in the quantum theory of the Coulomb field; (ii) determining the eigenstates of 
I the energy operator and assigning energy to the excitations of the theory; and (iii) proposing a 

, simple theoretical scheme to estimate the effect of the quantum fluctuations of the Coulomb field 

, ^ ■ on the energy levels of hydrogen-like atoms. We argue that the recent experimental advances in 

^ I hydrogen and muonic-hydrogen spectroscopy may provide the unique window of opportunity for the 

, verification of the Staruszkiewicz's theory. 

^ ■ PACS numbers; 04.62.-|-v 
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I. INTRODUCTION 



The theory of the quantum Coulomb field (QCF) has been proposed by Staruszkievificz [1-3] (see also Ref. [4] for the 
O ■ discussion of this theory from a different angle and Refs. [5-7] for other approaches to the quantization of long-range 
"j^ [ fields). It supplements standard quantum electrodynamics (QED) by describing long-range quantum fluctuations of 
^ . the Coulomb field [8]. The principal need for the development of this theory is that even though the electric charge 
^ ' of every particle type could be arbitrary in the QED framework, all electric charges are multiples of a single quantum 
' S ■ with a stunning accuracy [9] . 
^i^' The studies of the quantum Coulomb field have been centered on charge quantization and universality as well as 
f~i \ on the extensive search for a distinguished value of the fine structure constant (see e.g. Ref. [10]). While the former 
efforts have succeeded, the Staruszkiewicz's theory predicts that charged particles carry an integer multiple of the 
I same elementary charge, the latter efforts are still ongoing. Indeed, even though it was found that the mathematically 
^vq . distinguished interval of the allowed values of a is < a < tt, the hunt for the theoretical determination of the exact 
J> ■ value of the fine structure constant a = 1/137.036 ... is unfinished. We remark that to the best of our knowledge the 
. Staruszkiewicz's theory provides the first theoretical framework in which such a hunt is at all possible. 
■ Our goal is to extend the QCF theory to study its experimentally-relevant implications. We discuss the electric and 
magnetic field operators of the QCF, define a set of distinguished states of the QCF that might be carried by charged 
particles, and develop a simple theory providing a basis for the quantitative studies of the shifts of the energy levels 
OO [ of atoms resulting from the quantum fiuctuations of the Coulomb field. The possibility that there may be additional 
. level shifts absent in standard QED is strongly suggested by the puzzling recent precision spectroscopy experiment 
on muonic hydrogen [11]. 



II. BASICS OF THE QUANTUM COULOMB FIELD THEORY 



b. . . . . . 

Cu , In this section, we will briefly summarize the QCF theory of Staruszkiewicz. The central object of the theory is the 
phase field S{x). It may be derived on the classical level in two ways. 

Such a degree of freedom is always present for charged fields. This is seen by writing the action for the electromag- 
netic field and its source in the form 



d^xF^^'^Ff,, + / d''xC[eAf, + d^,S, . . .], (1) 



1 

~167r 

where the phase S of the charged matter field appears in the gauge-invariant combination eA^ + 9^5', while the dots 
stand for other degrees of freedom of the charged matter field, e.g., its real amplitude, relative phases between the 
different components of a spinor, etc. 

To illustrate this point, one can consider the simplest example of a field theoretical system with the elec- 
tric charge: the Klein-Gordon theory discussed in this context in Ref. [12]. Its gauge-invariant action reads 
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-ikl d'^xFt^'F,,^ + i / d'^x\{di^ + ieAi^)(j)\'^ - m'^\(j)\'^. Writing the Klein-Gordon field (f) as i?exp(i5) one finds that 
this action equals 

~ 1^ / '^^^^'"'^ ^1^- [di'Rd^R + R^id^S + eA^){d''S + eA^") - m'^R^] , 

in agreement with Eq. (1). In this particular case, the dots in the argument of C stand for R (the real amplitude of 
the charged matter field) and its derivatives. 

Coming back to the discussion of Eq. (1), we note that its variation with respect to A^ defines the current to 
be equal to -^d'^F^i,. The variation of the action with respect to S leads to the charge conservation law: d^j^ = 0. 
The momentum canonically conjugated with the phase field S is 

- _ _h 

The canonical quantization [S{x),Tts{y)]x°=y° = iS{x — y) results in 

[jo{x),S{y)]^o^yO = ie5{x-y). 

Integrating it over the hyperplane x^ = 2/°, one obtains 

[Q,S{y'>,y)]=ie, (2) 

whore Q = J d'^x jo is the charge operator, e = l/\/137.036 . . . stands for the unit of the electric charge, and h = c = 1. 
This interesting relation, however, is rather useless unless further input is provided about the phase field S{x). This 
is done in the following way. 

The electric charge can be also defined by looking only at the electromagnetic field. In order to find the minimal 
framework for such a consideration, one investigates the following scalar functional of the electromagnetic potential 
A^,{x): 

s{x) = -ex^'Ai^{x). 

It can be shown that it characterizes the electromagnetic fields of the Coulomb type completely. Namely, all electric 
and magnetic fields falling off as l/r^ can be uniquely cxprcsscid as thc^ appropriate derivatives of s(.t) through the 
relation —ex^F^,, = d^5{x) (see e.g. Ref. [9]). Moreover, s(x) satisfies the d'Alembert equation: Os = 0. 

In Ref. [1], Staruszkiewicz formulated a complete field-theoretical system by refiecting on the Gauss law. He 
identified S{x) (the property of the charged matter) with s{x) (the property of the electromagnetic field), and used 
the commutation relation (2) as the basis of the quantum field theory of the Coulomb field. Note that the identification 
of a degree of freedom of matter with some degree of freedom of the electromagnetic field is necessary if matter is 
required to carry a true electric charge even in asymptotic future/past (see e.g. Ref. [7]). 

The quantization of the phase field proceeds in two steps. Firstly, as the consequence of both the d'Alembert 
equation satisfied by S{x) = 3{x) and the required fall-off condition for S{x), one obtains the appropriate field 
operator [1]: 



S{x) = So - eQta,nh.{tp) + ^ ^ [c£mf£m + h.c] , 



^=1 TO=- 



2^1 (-^, i ; I ;tanhV) ^ - ^ tanh(^)2i^i (-f , ^; |;tanhV) 



where (r,0,(p) are the spherical coordinates around the charge, = arctanh(t/r) and 



^ + ir(^)- 

The phase field S is quantized outside of the light cone x^x'^ = t"^ — x"^ =0, i.e., for r > \t\. The non- vanishing 
commutators are 



[cim, c\,^,] = ATTe^5u>5mm', [<9, So] = ie. 
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As the second step of quantization, the representation space for the theory is defined to be the Fock space constructed 
upon the "vacuum" state |0), which is supposed to be annihilated by all C£m and Q (the last relation being highly 
non-trivial; see also [13]). 

The "simplest" state carrying n units of the electric charge, |n), is created from vacuum via exp(— m^o): 

|n) = exp(-m6'o)|0), Q\n) = ne\n), {n\m) = Snm- (3) 

We will refer to |n) as the n-charged vacuum state as it is annihilated by all C£„i. One should remember that |n) 
is a short cut for \n;{nem}), where ncm, the occupation of the {i,m) modes, is set to zero. The {i,m) modes are 
responsible for the global angular distortions of the Coulomb field. Another distinguished charged state of the QCF 
was found in Ref. [14]. 

Finally, we would like to draw attention of the reader to the fact that independently of the beautiful physical 
relevance attached to the field S{x), the quantum theory of this field is a perfectly well-defined, highly interesting, 
and well worked-out quantum field theory in 2 + 1 dimensional de Sitter spacctime (see e.g. the recent non-trivial 
results on the structure of the boost operators and the spectral decomposition of the vacuum state [15, 16]). 



III. FLUCTUATIONS OF THE QUANTUM COULOMB FIELD 



It is convenient for our calculations to switch to the following notation. For m = 0, we define 

Ceo = V^ie beo, Yeo{6, tp) = Y(o{d, ip), 



while for m > 0, we define 



27rie bfm ± b. 



'l,-m \ , 



Yem{e,<p)±Ye,-m{e,<p) 
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The redefined operators satisfy [^^m^^^on'] ~ fe'^mm/ ^'^nd , ^ 



= 0. 



From a given S{x) = s{x) the corresponding electromagnetic Coulomb fields can be computed using 

1 ~ 

A^'ix) = S(x). 

exx 

In the standard spherical orthonormal tetrad (N, 0, *): 

N = (sin 6 cos y , sin ^ sin ip, cos 9) , 
= (cos 6 cos ip, cos 9 sin ip, — sin 9), 
* = (— sin cos 0) , 



we find 



E = 



4 



t 



H= - 



1 1 



d 



t2 \sm9d(p " 89 



89 



S. 



Within the theory of Staruszkiewicz, the Coulomb electromagnetic fields are given only in the spacetime region 

\t\ < r (the theory is originally formulated in spatial infinity outside of the light cone, where this condition does 
not cause problems; we extrapolate the fluctuations from the spatial infinity to near the charge center assuming 
homogeneity of the Coulomb field). This refiects the lack of a complete theory of charged matter accompanied 
by such fields (for attempts in this direction see Ref. [7]). In this paper, we circumvent that difficulty by dealing 
exclusively with effects on the surface t = 0, for which the theory gives complete information about the electromagnetic 
fields (apart from the singular point r = 0): 



E(0,x) = 



H(0,x) = 




-hh.c. 



N 



_ ^ 

_sm9 8ip 89 



(4) 
(5) 



e=i m=-e 
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Several remarks are in order now. 

First, these expressions show that there is a fluctuating electric and magnetic field around a charge: an expected 
feature of the quantum Coulomb field. The amount of fluctuation depends on the quantum state of its Coulomb field. 
Physically sound states describing the QCF of a charge at rest are those reproducing the classical result: 

(E(0,x))=ne^, (H(0,x))=0. (6) 

This condition, however, is satisfied by the infinite number of states. For example, by any single Fock state \n, {n^m})) 
where nem is an arbitrary non-negative integer. It is therefore fundamentally important to ask which quantum state 
of the Coulomb field is carried by a charged particle, say an electron or a proton ? 

Second, the amount of the fluctuation of the quantum Coulomb field should be experimentally measurable through 
the studies of the level shifts of atoms (especially hydrogen-like for which the techniques of precision spectroscopy 
work best [17, 18]; we will come back to this point later). This is a simple conclusion well founded on the success of 
QED to explain the Lamb shift. The Lamb shift, originally the splitting of the 2Si/2 and 2Pi/2 levels in hydrogen, 
results from the coupling between the orbiting charge and the fluctuating vacuum fleld [19-21]. 

Third, one easily sees from (4) that the relative fluctuation of the QCF with respect to its classical counterpart 
is (i) independent of the distance from the charge center as both decay as 1/r^ and (ii) decays as the inverse of the 
charge carried by the QCF (fluctuations arc charge independent while the classical fleld goes linearly with charge). 
The (i) effect is in stark contrast to what one finds in QED, where the relative modification of the classical Coulomb 
field by virtual electron-positron pairs is short-ranged [21, 22]. 



IV. ENERGY ASSOCIATED WITH THE QUANTUM COULOMB FIELD 

Having realized that there are many quantum states of the Coulomb field that are physically allowed, we propose 
a simple physically-motivated approach allowing for finding a set of distinguished states of the QCF. Namely, we 
propose to "label" different configurations of the QCF by an energy associated with them. To this aim, we define 
electromagnetic energy operator on the surface t = 0: 



£c 



8nJ 



E(0, x)E(0, x) + H(0, x)H(0, x) :, 



(7) 



where :: stands for normal ordering [23]. 

In classical electrodynamics the expression (7) corresponds to the energy of field configurations. In the QCF theory 
states reproducing on average the classical Coulomb field do exist, and we feel it justified to claim that their energy 
should be computed from the quantum version of the classical energy functional. However, since E(0,x) and H(0,x) 
scale as 1/r^, this operator has problems due to the divergence of the integral at r = 0. Moreover, it does not have 
the interpretation of the generator of the time-translation symmetry (the Hamiltonian). We will discuss the former 
problem below and mention here that it is our assumption to assign the energy to the excitations of the QCF through 
the operator (7). 

Putting operators (4) and (5) into the above expression, we obtain 



e=i 
A(-l) 



We introduced above 



{bl + bll) + E ^IJ^rn + bl.Je,-m + (-)'"+^f (^L " ^l-^ + h-c.) 



m=l 



4 - Gji{£ + 1) _ 2_ Gltji + 1) 



(8) 
(9) 



and 



f (fix 

Stt 



/ 



(10) 



The latter corresponds to the classical electromagnetic energy of a unit charge at rest. We do not propose any new 
solution to the obvious lack of convergence of this integral at r = 0. On general grounds, however, we expect that 
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in any such solution the parameter Sc should be of the order of the electron's rest mass m,.. This is in line with the 
century long efforts to associate the mass of an electron with the energy of its field [24] (see also Ref. [25] for the 
relevant hadronic example) . We shall also argue below that the actual spectroscopic consequences of the quantum 
theory of Coulomb field will be in a good position to determine the numerical value of Ec- 

To look for the spectrum of the electromagnetic energy operator, we diagonalize (8). This is done by applying the 
Bogolubov transformation 



where assuming that m > 



Bim = cosh{ui,n) 'bem " sinh('u<!m)S! , (11) 



(-)^+iarctanh(/3^) T(-)™+^arctanh(/3^) 
um = ^ ,ue,±m = ^ • (12) 



The transformation (11) can be easily expressed as the explicit function of £ employing 

2 

arctanh(/3£) = In 



Finally, we note that [Bem, bI,^,] = Su/S„ijn' and [B^^, -6^'^/] = hold for the redefined operators. 
Putting Eqs. (11) and (12) into Eqs. (8) and (9), we find 



where we have employed the identity — = 2^£{£ + 1) to simplify the result. 

The first term in (13) recovers the classical electromagnetic energy of a charge at rest. The second one is the 
contribution of the zero point modes (see also [23]) 

S E 7. - l) {2i + 1) « -5 X 10- J. (14) 

The third term in Eq. (13) provides the spectrum of the electromagnetic energy operator, which is degenerate for all 
m in every £ sector. We define the ground state of £c, say \G), as the state that is annihilated by all Btm operators. 
For generality, we assume that it carries n quanta of the electric charge: Q\G) = ne\G). The excited states are 
created in the n-charged sector by acting (an arbitrary number of times) Bj^ on It is worth to mention that the 
excitation gap to the £ sector equals 

2ecV£(£ + T) 

which is a huge number assuming that indeed Sc = 0{me). Given the physical interpretation of the electromagnetic 
energy operator and an expected huge gap in the excitation spectrum, we propose the ground state \G) as the natural 
candidate for representing the quantum state of the Coulomb field of a stable charged particle (e.g. an electron or a 
proton [26]). Note that both the charged vacuum state |n) and the ground state |G) are spherically symmetric. The 
expectation value of the energy operator is lowest in the state \G) (the state jn) is not an eigenstate of this operator). 
In what follows, we describe some of its basic properties. 

First, we note that the expectation value of the electric and magnetic field operators, calculated in the ground state 
[G), reproduces the classical result (6). States supposed to correspond to moving particles can be constructed because 
the generators of the boosts are explicitly known (see Refs. [27, 28] and the Appendix). For example, the quantum 
state of the Coulomb field of the charge moving in the +z direction with velocity v =tanhA is 



|G',A) =exp(-iAM"^)|G) 



where is the boost operator (A6). 
We expect that the averages. 



(G,AlE(0,x)lG,A), (G,AlH(0,x)|G,A), 
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computed as the functions of x on the surface t = 0, are equal to the classical fields of a moving charge crossing this 
surface at x = [29] : 

ne 1 — 

H (1 - ^2 sm"^ 9)^/^ 

The theory has enough structure to prove that it is the case, but as a check we have verified the result by a direct 
calculation to the order A^. We also note that the action of the boost operator on the ground state \G) additionally 
populates the {£, m) modes in the jG, A) state. This allows for the addition of the non spherically-symmetric component 
to the quantum Coulomb field, which is indispensable for the reproduction of the above classical result. 

Second, it turns out that \G) is a squeezed vacuum state \n) (see Ref. [30] for the discussion of the squeezed states 
in the quantum optics context): 

oo t 

\G) = exp(D)|n), ^ = 2 E E (^J™ " ^'m) • (16) 

^=1 m=-e. 

This is easily seen from the following identity: exp(— D)^^^ exp(L') = b^m- Moreover, the overlap between the ground 
state and the charged vacuum state is 



{n\G) = n 



cosh 



/ arctanh(/3£) 
V 2 



0.997. 



The two spherically-symmetric states are therefore quite similar. Interestingly, their overlap is independent of the fine 
structure constant a = e^. Thus, it is of purely geometric nature. Furthermore, we can study the population of the 
(^, m) modes, 

It equals about 3 x 10""^ for i = 1, decays motonically with (., and for £ 1 it approaches 1/64^^: only low {l,m) 
modes are noticeably populated. 

Third, while the expectation value of the field operators (4) and (5) reproduces the classical result, their fluctuations 
in the ground state have absolutely no classical counterpart. This is so because the long-range fields are now described 
by a quantum field theory, and not just by the classical functions. Indeed, 

(E(0,r,e,(p)E(0,r, 6*', </?')) = — ^ + -7> 'niPeicosecose' + cosUp - ip')sm.6s.m.e'), 

f f . Gf 

£=1 ^ 

where % = 1 when the average is calculated in the n-charged vacuum state |n) and 



when it is calculated in the ground state {Pi stands for the Legendre polynomial.) This shows explicitly that 
the charged states |n) and \G) (i) differ by the pattern of the fluctuations of the electric fleld associated with them; 
(ii) are spherically symmetric. 



V. DISCUSSION OF THE RESULTS 



To start, this work proposes that there can be a fluctuating long-range electromagnetic field around every charge. 
This field originates from the quantum fluctuations of the Coulomb fleld. We propose that it should be considered 
in addition to the fluctuating vacuum field studied in the standard QED framework. Since such a field has not been 
observed yet, the confirmation of its presence will fundamentally update our understanding of charged particles. It is 
thus critical to find out the experimentally accessible consequences of this conjecture. 

General experience with quantum systems interacting with the fluctuating electromagnetic fields is that (at the 
very least) their energy levels are subjected to shifts [21]. We propose to employ the spectroscopic measurements 
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of the energy levels of the bound systems such as hydrogen, muonic hydrogen, muonium, etc. in the hunt for the 
experimental evidence of the quantum nature of the Coulomb field. 

To quantify the expected level shifts, one can adopt the perturbative formalism analogous to the one employed 
for the calculation of the Lamb shift [19-21, 31]. In this approach the quantum state of the Coulomb field of a 
nucleus is modeled by our ground state \G) (16). The field operators used for writing down the coupling between 
the orbiting charge and the fluctuating Coulomb ficild are given by Eqs. (4) and (5). The excited states of the QCF 
and their energies (both required in the perturbative expansion) are provided by the eigenstates and eigenvalues of 
our electromagnetic energy operator (13). Since the eigenvalues depend on Sc, we propose to treat it as the (only) 
free parameter and argue that its value can be experimentally determined from the spectroscopic measurements. 
Naturally, the larger the smaller the level shift due to the QCF should be. 

The magnitude of the level shifts resulting from the presence of the QCF is bounded by the discrepancies between 
the experimental measurements and the QED calculations. For example, the frequency of the 1S-2S transition in 
hydrogen is experimentally known with the fantastic accuracy of about 40 Hz corresponding to the relative accuracy 
of about 1 part in 10'^^ [17, 18]. Even more amazingly, this result is expected to be significantly improved in the 
foreseeable future [18]. The QED calculations match it with the relatively thick error bar of several tens of kHz [22]. 
The most of this theoretical uncertainty comes from the imprecise knowledge about the "distribution of charge" in a 
proton. Thus, the QCF corrections to this transition are bounded from above by a few tens of kHz. 

The above-outlined possibility to consider additional level shifts in hydrogen-like atoms is especially appealing in 
the light of the current discrepancies between the standard QED predictions and the spectacular recent measurements 
of the Lamb shift in muonic hydrogen [11]. It is also quite natural given the fact that precision spectroscopy has 
been serving for a long time as the Rosetta Stone for deciphering the laws of quantum physics [18]. The formalism 
proposed in this manuscript should lay the ground for the studies of the spectroscopic consequences of the quantum 
Coulomb field. The work along these lines is already ongoing [31]. 
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Appendix 

We use spherical harmonics defined as 



Yem{0,y^) = (-)(-+l-l)/V^ ^^^+ ^ Hj; pl-l(cosg)exp(im^), (Al) 

following the convention introduced by Staruszkiewicz in Ref. [3]. For the convenience of the reader we list the first 
three £ sectors: 

1 / 3 

Yoo = -7= , Yio=i\ — cos 9, 
V47r V 47r 



Yi,±i = Ti\J^ sine e^"^, = y/_i-(i _ 3 cos^ 6), 



Y2,±i = ±;/^ sin(2^)e±'^ F2,±2 = -]/ ^ sin^ ee^''^. 



In particular, this implies that 



YeUO, = (-)"+'r^,-„^(0, L±Y,^ = V^(^+l)-m(m±l)r,,™±i, (A2) 

where 

In our manuscript, we use the redefined spherical harmonics Yim- Considering m > 0, one has 

YL{9,^) = i-r+'YUe,^), Yl.mif),^) = -{-r+%-m{0,^). (A3) 

Noting that 



^ Yem{0,^)Y;^{9',^')= Ycm{9,^)Yl^{9',^'), 



m=—i m=—t 

we obtain 

t 

^ Yo^(9.in\Y* (9'.in'\ = - 

4n 



V Yem{0,y^)Y;(9',ip') = ^^-tlp^ (cos 61 cos 61' + cos{ip - ip')sm9sme'), (A4) 

m=-i 
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where Pi is the Legendre polynomial. Eqs. (A2), (A3) and (A4) are useful in the derivation of both the electromagnetic 
energy operator and the expression for the two point correlation function. 

The boost operator M°^, that we use in the paper, was calculated in Ref. [27] (see also Ref. [28]). We list it below 
for reader's convenience: 

= TE^e^^^'' + + 4^^ E E ^ ^p^^y^^^i4mce-i,m - h.c.), (A5) 
which in our notation corresponds to 

= ^\/J^(^io - bl) + i £ E \/5^v/^^^(^L^^-i,m - 11.C.). (A6) 

We mention in passing that the derivation of the boost operators (A5) and (A6) assumes that spherical harmonics 
are given by Eq. (Al). 



